Two infinite families of nonadditive quantum error-correcting codes 
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We construct explicitly two infinite families of genuine nonadditive 1 -error correcting quantum codes and 
prove that their coding subspaces are 50% larger than those of the optimal stabilizer codes of the same param- 
eters via the linear programming bound. All these nonadditive codes can be characterized by a stabilizer-like 
structure and thus their encoding circuits can be designed in a straightforward manner 



One major family of quantum error-correcting codes 
(QECCs) III El El ID, which are powerful tools to fight the 
quantum noises in various quantum informational processes, 
are called as additive or stabilizer codes Q |6l |2l jH- The 
coding subspace of a stabilizer code is specified by the joint 
+ 1 eigenspace a group of commuting multilocal (direct prod- 
uct of) Pauli operators. Usually [[n, k, c?]] denotes a stabilizer 
code of length n, the number of physical qubits, and distance 
d, i.e., correcting up to J -qubit eiTors, that encodes k log- 
ical qubits (2'^-dimensional subspace). 

The first example of nonadditive codes, codes that cannot 
be described within the framework of stabilizer, was an infi- 
nite family of 1 -error-detecting codes f9','T0|, e.g., ((5, 6, 2)), 
with coding subspaces being 50% larger than the optimal sta- 
bilizer codes of the same parameters. Recently another family 
of 1 -error detecting codes with still larger encoding subspaces 
has been constructed in [11, 1 and slightly improved in |12|. 
Here we have denoted by ((n, K, d)) a nonadditive code of 
length n and distance d that encodes a if-dimensional logical 
subspace (about log2 K logical qubits). 

The first example of nonadditive code |fT3ll . namely 
((9,12,3)), that outperforms all the stabilizer codes of the 
same length while capable of coiTecting arbitrary single qubit 
errors has recently been constructed via a graphical approach 
based on graph states. Later on an optimal 10-qubit code 
((10,24,3)) has been found via a comprehensive computer 
search lfT2l . Recently a family of codes of distance 8 that 
encode 3 more logical qubits than the best known stabilizer 
codes have been constructed based on nonlinear classical 
codes fl4i. However the possibility of being equivalent to 
some subcode of an optimal stabilizer code or even a stabilizer 
code of the same parameters under local unitary transforma- 
tions has not yet been excluded. 

Generally, being without a stabilizer structure, the nonad- 
ditive codes promise larger coding subspaces while they are 
harder to construct and identify than the stabilizer codes. On 
one hand there is no systematic construction so far and all the 
good codes are found via computer search fT2l[T5ll , which is 
impossible for a relatively large length (e.g. n > 11). On the 
other hand an obvious criterion for a genuine nonadditive code 
is to check whether or not its coding subspace is larger than 
all the stabilizer codes of the same parameters. However the 
exact bound for stabilizer codes is generally unknown. As a 
result it is of interest to find nonadditive error-correcting code 



when the length tends to infinity that outperforms all the sta- 
bilizer codes of the same parameters. 

In this Letter we shall construct two infinite families of gen- 
uine nonadditive 1 -error-correcting codes with coding sub- 
spaces being 50% larger then the corresponding optimal 1- 
error-correcting stabilizer codes of the same parameters to 
show that the nonadditive error-correcting codes outperform 
the stabilizer codes even when the length n tends to infinity. 
All the nonadditive codes are characteiized by a stabilizer-like 
structure and therefore the encoding-decoding circuits can be 
designed in a straightforward manner. 

Two families of nonadditive 1 -error coiTecting codes that 
we shall construct have the following parameters 

D(m,a) = ((K.,i2^""'""',3)) (1) 

where = ^'""g"^''^ + a with a = 0, 1 and m > 1. To 
ensure that they are genuine nonadditive we shall prove that 
the corresponding optimal 1 -error-correcting stabilizer codes 
of the same length have parameters 

[[7V;^,7V^-2™-6,3]] (2) 

by working out analytically the linear programming bound for 
the lengths N!^. Notice that the quantum Hamming bound 
permits exactly one more logical qubit, i.e., [log2(3iV^ + 
1)] = 2m + 5. The first nonadditive codes of these two fami- 
lies are ((41, 3-232, 3)) ((42, 3-233, 3)) respectively while 
the optimal stabilizer codes have the parameters [[41, 33, 3]] 
and [[42,34,3]]. 

Our construction is based on a family of stabilizer codes of 
lengths {22''+3}jj^^ jjllsl and two nonadditive codes of length 
9 and 10 discovered recently lfT2l|T3l and is a kind of pasting 
stabilizer codes with nonadditive codes that generalizes the 
pasting of stabilizer codes in Ref.fT6l. We denote by Xy^ 
three Pauli operators acting nontrivially only on some qubit 
labeled by v and by I the identity operator. Furthermore for 
a given index set U we denote Xjj — Htjec/ ™^ similarly 
for other Pauli operators. 

Let us look at the optimal stabilizer of length 2'^'^^^ at first. 
According to Ref.|i5J the stabilizer of the code has 2r + 5 gen- 
erators with two of them being Xu^ and Zu^ where we have 
labeled 2'^''+^ physical qubits with f/^ = {1, 2, 3, . . . , 2"^-"+^}. 
The remaining 2r + 3 generators are given by 

[SI = ;fh.2h'=-i+hi+h2,.+3 \ kC,Ur]. (3) 
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Here denotes a 2^''+^-dim vector that is the A:-th row 
of a (2r + 3) X 2^*"+^ matrix i/,. = [cq, ci, . . . , C22r+3_i] 
whose fc-th column Ck being the binary representation of k, 
e.g., = (0,0,..., 1) and C22r+3_i = (1, 1, . . . , 1) and 
ho = is the zero vector And for a vector h with com- 
ponents {hy I V E Ur} we have denoted = Yiveu '^v" 
andZh = n,^^^4'^. 

Despite their nonadditiveness the codes ((9,12,3)) and 
(10, 24, 3)) admit a stabiHzer-Hke structure and can be most 
conveniently formulated by using the graph states ||20l \2T\ . 
We denote by G = {V, E) a simple undirected graph with a 
set V of vertices and a set E of edges. Two vertices are con- 
nected with an edge iff {a, 6} e E. Two graphs Ga (a = 0, 1) 
on |Vb| = 9 and \Vi\ = 10 vertices are shown in Fig.l. By 
labeling \V\ qubits by V we can define the graph state corre- 
sponding to the graph G a |G) = Ug\+)'^ where 



ua= n 

{a,b}eE 



1 + Za + Zb — ZaZb 



(4) 



and denotes the joint +1 eigenstate of Xy for v G V. 
Obviously Uq = 1 and the graph state |G) is also the +1 joint 
eigenstate of the following n stabilizers = UgXJJg- 

From the graph state |G) a basis for the whole system a 
basis of the whole system can built {Zc\G) \ C C V}. For 
trivial graph with no edge the graph-state basis reads Zc\+)X ■ 
Any given collection of vertex subsets of V will define a basis 
that spans a subspace which is referred to as the graph-state 
basis. With a graph and a collection of vertex subsets will 
define a subspace. 



For the nonadditive code ((9, 12, 3)) 



"(0,0) 



we consider 



the loop graph Go = (Vb, ^'o) on 9 vertices as shown in Fig. la 
and corresponding graph state |Go). The joint +\ eigenspace 
of the following 6 observables 



ai 


= ^Go'^{3,8}^Go7 


(5a) 


0-2 


= ^Go'^{6,2}^Go: 


(5b) 


az 


= ^Go'^{9,5}^Go7 


(5c) 


Ax 


= ^Go ('^{4,7,3,6,9} V69)WGo J 


(5d) 


A2 


= ^Go('^{l, 7,3,6} "^39)^00 ) 


(5e) 




= ^Go ('^{1,4,3,9} V36)^Go J 


(5f) 



is exactly the 12-dim coding subspace of the code D(o_o) ■ Here 





FIG. 1: Two graphs Ga 
{1,2,..., 9} and Vi = {0,1, 



{Va,Ea) are shown with Vo = 
. , 9}. They provide the graph-state 



we have denoted Vab = \{^ + Xa + X^ — XaXf,). If we denote 
^0 = -41^2-^3 we have Try^^o = 2^. It should be noted 
that the code is nondegenerate, which can be easily seen from 
its weight distributions. 



The nonadditive optimal code ((10, 24, 3)) = 



''(0,1; 



has a 



graph-state basis corresponding to the graph Gi — (Vi, Ei) 
on 10 vertices as shown in Fig. lb and the corresponding graph 
state is |Gi). Obviously the graph is unchanged under the 
following two permutations 

1T = (14)(23)(69)(78), r = (12)(34)(67)(89). (6) 

that act naturally on V . Two swapping operators A^^r and 



M.r are defined via, e.g., M.-nZc\ 



fori 



arbitrary C QV . For later use we define two controlled-swap 
operations with qubits and 5 as sources respectively as 



(1 



1 



Xq 



(1 
(1 



A'o)XO' 



(7a) 
(7b) 



which can be readily implemented via Toffli and control-not 
gates. According to Ref.fTZ). those 24 graph-state bases of 
the code D(o,i) are given by {|G^j^) = Z^i |Gi)} with i = 
1, 2, . . . , 6 and /i, ^ = 0, 1 where 24 subsets of Vi are 



T^id) AvBA tiT^iA) 



(8) 



with AAB = AU B — AnB being the symmetric difference 
between two sets A and B and 

A={O,2,3},B = {5,l,2},Gi = 0, 
G2 = (1, 2, 3, 9}, G3 = (1, 2, 7, 8}, G4 = (1, 2, 6, 7, 9}, (9) 
G5 = {1,3,7,8,9},G6 = {1,3,4,6,7,9}. 



It turns out that the nonadditive code 



''(0,1: 



admits also a 



stabilizer-like structure. With the help of the following en- 
coding operator 



^cnc — Z2UgiTt%;Z2 



(10) 



it can be readily checked that the following 6 mutually com- 
muting observables 



/3i 

Bo 
B2 



Wcnc'^{2,3,7}^cncJ 
^onc'^{6,7,8}^cnc' 
^onc'^{3,4,6,9}^enc) 
^cnc ('^6V67)Z^lnc: 
Wonc (-^{1,2} V37) Z^on 
^onc (A'4V36)Z^cnc 



(11a) 
(lib) 
(11c) 

(lid) 
(lie) 
(llf) 



stabilize the coding subspace of D(o,i), i-C-, one one hand 
Pk\Cl^) = Bi\Cl^) = \Cl^) for all possible k, I, i, fi, v and 
on the other hand the joint H-1 eigenspace of these 6 observ- 
ables, whose projector is given by 



bases for two nonadditive codes ((9, 12, 3)) and ((10, 24, 3)). 



B,. 



(12) 
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TABLE I: The stabilizing observables of thie nonadditive codes DJ^ 
whose piiysical qubits are labeled with (7m U. . .U?7iUVa (a = 0, 1) 
The blank entries stand for suitable identity operators lu^. or ■ 







Um-l 


■ (72 Ui VoorVi 


0(0,1) 
0(0,1) 




0(0,1) 
0(0,1) 


sr 

om 
'-'2 






0(0,1) 


c-m 


^m-l 






Om 

'->2m-6 


— 1 
'-*2m-S ' 




^0(0.1) 
Ai{0,l) 


cm 

'-'2m-5 
om 

'->2m-4 


om — 1 
'-'2m-7 ■ 

om— 1 
'-'2m-6 ■ 


• 



^2m-l 

^o(o,i) 



5 m 
2m-3 

5 m 
2 m -2 



■m — 1 
2m — 5 

m— 1 
2m-4 



5i Xui 



51 ^C; 



>-^2m + l 

'^2m + 2 

0(0,1) 



^2m + 3 
(0,1) 
2m+4 

*-^2m + 5 



5 m 
2m-l 

om 

5 m 
2m + 2 

om 

<->2m+3 



m— 1 
2m -3 
m— 1 
2m — 2 
m— 1 
2m— 1 

;7m — 1 

''2m 



om — 1 
*-'2m + l 



Si 

s! 



^6 

5? 



51 
51 



Qfi or /3i 
"2 or /32 
a-i or /33 
^1 or Bi 
A2 or Z32 



(0,1) 
2m + 6 



Ao or So 



has exactly dimension 24, i.e., TrviV = 24 since Trv^/Bo — 
2^. An encoding circuit can therefore be designed in a simi- 
lar manner as that of ((9, 12, 3)) 1 13J . We note also that this 
nonadditive code is non-degenerate. 

Now we are ready to present our construction. We consider 
qubits and label them by disjoint set Um U Um-i U . . . U 
Ui U Va with \Uk\ = 2^'"+^ and |Vq| = 9 + a with fc < m and 
a = 0, 1. We claim that the joint +1 eigenspace of those 2m + 
6 observables {Ol"^}^^^^ with a or 1 as defined in Table 
I is the code D(„j a) in Eq.([Tj) with the following projector onto 



the coding subspace 



2m+6 



-pa ^ 



n - 



o 



(a) 



(13) 



In Table I observables 5f are defined in Eq.pl and {a^, Aj}, 
{Pi,Bj} are defined via Eq.ljs]) and Eq.( 11 1 respectively. 
Blank entries represent suitable identity operators. By jux- 
taposition of some operators in the same row we mean their 
direct product. 

First of all, these 2m + 6 stabilizing observables detect all 
2-qubit errors because firstly all errors happened on U -blocks 
or V blocks can be detected because all the subcodes are pure 
1 -error correcting codes and secondly all two errors happened 
on different qubit blocks can always be detected by the sta- 
bilizer containing Xjj^ and Zu^. for some k. Thus we obtain 
a pure 1 -error correcting codes of length N^. Secondly, by 



noticing TrO 
we have 



(a) 



Ofori < 2TO + 5andTra 



(a) 

2m+6 



Tr(l + 0^1+6 

22m+6 



^ 2 



(14) 



Thus we obtain the 1 -error correcting code of parameters ex- 
actly as given in Eq.Q. Now we shall demonstrate that its 
coding subspace is 50% larger than the corresponding opti- 
mal stabilizer codes so that our codes are genuine nonadditive 
codes that are neither equivalent to some stabilizer codes un- 
der local unitary transformations nor subcodes of some larger 
1 -error correcting stabilizer codes of the same length. 

The quantum Hamming bound for a 1 -error correcting sta- 
bilizer code, e.g., n—fc > [log2(3n+ 1)] for a stabilizer code 
[[n, fc.3]], being introduced initially for the non-degenerate 
codes, is valid for both degenerate and generate codes of dis- 
tance 3 and 5 i6J and of a large enough length |fT9]| . In the 
case of rt = we have the quantum Hamming bound 
n — k > 2m + 5. This is not enough to prove the nonad- 
ditiveness of our codes. However by working out analytically 
the linear programming bound we have 

Theorem If there exists a stabiUzer code [[A^^, k, 3]], de- 
generate or non-degenerate, with iV^ = - — g — - + a and 
m > and a = 0, 1 then iV,° - A: > 2m + 6. 

Proof. Given a stabilizer code [[n, fc, d]] its weight distribu- 
tions {Ai}"^Q is defined by 



A,. 



1 

22fe 



{i = 0, 



(15) 



where the summation is over all errors supported on i qubits 
and P is the projector onto the coding subspace. It is obvious 
that A^ > 0, Ao 1, and At = 2^* so that {Aj/2'} can be 
regarded as a probability distribution with s = n — k. For an 
arbitrary function f{x) we denote its average 



(16) 



In the following we shall formulate a subset of the linear pro- 
gramming bound for 1 -error correcting code, which serves our 
purpose perfectly. For a complete set of linear programming 
bound see Ref.llglfTTl. 

Linear Programming bound (Restricted) If there exists a 
stabilizer code [[n, fc, 3]] then the following conditions hold 
true 



Ai = (3?i - 4a;), 
^((4a;-3n+l)2-3n-l), 
LtJ 

£^2, >2^-i. 



(17a) 
(17b) 

(17c) 



In the case of a = 0, i.e., n = with m > we introduce 
a nonnegative function f{x) = (3n + 1 — Ax)^ and it is easy 
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to check that as long as n > 5 

/(O) = (3n + l)2 > (3n + 5)(3n-7) + 16, (18a) 

/(I) = (3n-3)2 > 4(3n + 5), (18b) 

/(2) = (3n-7)2 >2(3n + 5) + 16. (18c) 



If there exists a stabilizer code [[n, k, 3]] then Eqs.( 17a|17c i 
must hold. As a result 

(/(x)) =3^ + 1 + 4^1 + 2^2, (19a) 
LfJ LtJ 
16 + 16^2 + ^ /(2i)yl2» > 16^^2^ > 8.2", (19b) 



i=2 



i=0 



where we have used that fact that /(2i) > 16 since '^^^^ , the 
unique zero of f{x), is an odd integer Putting all these pieces 
together we obtain 



1=0 



LfJ 



> 



/(O) + /(l)Ai + /(2)A2 + f{2i)A2 



i=2 



> /(O) - 16 + f{l)A, + (/(2) - 16)v42 + 8.2" 

> (3n + 5)(3n-7 + 4Ai+2A2) + 8.2" 

= (3n + 5)(/(x) - 8) +8.2", (20) 

in which the strict inequaUty comes from the /(O) term. Tak- 
ing into account of {f{x)) > 8 we obtain 2" > 3n + 5, i.e., 
n — k > 2m + 6. 

In the case of a = 1, i.e., n — N^^^ with m > we define 
g{x) = (3n + 2 — 4x)(3n — 2 — 4a;) which is nonnegative 
on integers because ^"^^ is an integer. It is obvious that as 
long as n > 5 we have g{i) > 2(3n + 2) for i = 1,2 and 
most importantly g{0) > (3n + 2)(3n — 4). If there exists a 
stabilizer code [[n, k, 3]] then Eqs.( 17a|l7c i must hold, which 
leads to {g{x)) = 3?! — 4 + 2Ai + 2A2. As a result we have 

2Hg{x)) > 5(0) + ff(l)Ai+ 5(2)^2 

> (3n + 2)(3n-4 + 2yli +2A2) 

= (3n + 2)(5(x)), (21) 

in which the strict inequality sign is due to the g{0) term. 
Since {g{x)) > we have 2" > 3n + 2, i.e., n - k > 2m + 6. 
□ 

It should be noted that the optimal stabilizer codes of pa- 
rameters as given in Eq.Q exist and construct is already given 
by the stabilizers in Table I with the stabilizers acting on 
qubits Vo or Vi being replaced by 6 stabilizers of the pure 
optimal stabilizer codes [[9, 3, 3]] or [[10, 4, 3]]. 

The stabilizer-like structures of our codes simplify signifi- 
cantly the encoding and decoding procedures. Let us suppose 
we have akeady the encoding and decoding circuits for the 
codes D(o.a) and for the Gottesman's codes. With some addi- 
tional controlled-not gates in front of the encoding circuits of 
these individual codes we obtain the encodings of our codes. 
To decode we have only to check at first the first 2m genera- 
tors in Table I, from which we can be sure wether the errors 



happen on some U-block or not. If yes we use the decodings 
for Gottesmans codes and if not then we have only to decod- 
ing D(o,o) with the detailed circuit in the case of a = being 
given in lfT2l . 

We acknowledge the financial support of NNSF of China 
(Grant No. 90303023, 10675107, and 10705025) and the 
A*STAR gi-ant R- 144-000- 189-305. 

Note added. On finishing the paper another infinite family 
of genuine nonadditive codes has been reported in I22J . 



[1] P. W. Shor, Phys. Rev. A, 2, 2493 (1995). 

[2] C. H. Bennett, D. P DiVincenzo, J. A. Smolin, andW. 

K.Wootters, Phys. Rev. A 54, 3824 (1996). 
[3] A. Steane, Phys. Rev. Lett. 77, 793 (1996). 
[4] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 
[5] D. Gottesman, Phys. Rev. A 54 1862 (1996). 
[6] D. Gottesman, arXive: quant-ph/9705052 
[7] A. Calderbank, E. Rains, P. Shor, and N. Sloane, Phys. Rev. 

Lett. 76, 405 (1997). 
[8] A. Calderbank, E. Rains, P. Shor, and N. Sloane, IEEE Trans. 

Inform. Theory, 44, 1369 (1998). 
[9] E.M. Rains, R. H. Hardin, PW. Shor, and N.J.A. Sloane, Phys. 

Rev. Lett. 79, 953 (1997). 
[10] E.M. Rains, IEEE Trans. Inf. Theory 45, 266 (1999). 
[11] I.A. Smolin, G. Smith and S. Wehner, Phys. Rev. Lett. 99, 

130505 (2007). 

[12] S. Yu, Q. Chen, and C.H. Oh, arXiv: 0709.1780vl [quant-ph] 
[13] S. Yu, Q. Chen, C.H. Lai, and C.H. Oh, Phys. Rev. Lett. 101, 
090501 (2008) 

[14] M. Grassl and M. Roetteler, Proc. 2008 IEEE Int. Symp. on Inf. 

Theory (ISIT 2008), 300 (Toronto, Canada, luly 2008). 
[15] A. Cross, G. Smith, I. Smolin, and B. Zeng, IEEE Trans. Inf. 

Theory 55, 433 (2009). 

[16] D. Gottesman, arXive: quant-ph/9607027 

[17] E. Rains, IEEE Trans. Inform. Theory 44, 1388 (1998); ibid, 

IEEE Trans. Inform. Theory 45, 2361 (1999). 
[18] P Shor and R. Laflamme, Phys. Rev. Lett. 78, 1600 (1997). 
[19] A. Ashikhmin and S. Litsyn, IEEE Trans. Inform. Theory 45, 

1206 (1999). 

[20] D. Schlingemann and R.R Werner, Phys. Rev. A 65, 012308 
(2001). 

[21] M. Hein, J. Eisert, and H.J. Briegel, Phys. Rev. A 69. 
062311(2004). 

[22] M. Grassl, P Shor, G. Smith, J. SmoUn, and B. Zeng, arXiv: 
0901.1319 [quant-ph]. 



